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ABSTRACT: 

We aimed to identify microRNA (miRNA) expression patterns in the serum 
of prostate cancer (CaP) patients that predict the risk of early treatment failure 
following radical prostatectomy (RP). Microarray and Q-RT-PCR analyses identified 
43 miRNAs as differentiating disease stages within 14 prostate cell lines and reflected 
publically available patient data. 34 of these miRNA were detectable in the serum 
of CaP patients. Association with time to biochemical progression was examined in 
a cohort of CaP patients following RP. A greater than two-fold increase in hazard of 
biochemical progression associated with altered expression of miR-103, miR-125b 
and miR-222 (p<.0008) in the serum of CaP patients. Prediction models based on 
penalized regression analyses showed that the levels of the miRNAs and PSA together 
were better at detecting false positives than models without miRNAs, for similar level 
of sensitivity. Analyses of publically available data revealed significant and reciprocal 
relationships between changes in CpG methylation and miRNA expression patterns 
suggesting a role for CpG methylation to regulate miRNA. Exploratory validation 
supported roles for miR-222 and miR-125b to predict progression risk in CaP. The 
current study established that expression patterns of serum-detectable miRNAs taken 
at the time of RP are prognostic for men who are at risk of experiencing subsequent 
early biochemical progression. These non-invasive approaches could be used to 
augment treatment decisions. 



1. INTRODUCTION 

In men in the USA and elsewhere, prostate cancer 
(CaP) is the most common noncutaneous cancer diagnosed 
and second leading cause of death [1, 2]. At multiple 



stages of disease there is ambiguity over progression risk 
and treatment responses. For example, the prostate-specific 
antigen (PSA) [3, 4] test has significantly increased 
cancer detection but has poor specificity and prognostic 
accuracy. The result of this ambiguity is that, amongst 
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men who have undergone initial therapy, it is uncertain 
who will relapse with recurrent aggressive disease [5, 6]. 
These ambiguities are clinically relevant because men who 
experience treatment failure have a significantly increased 
risk of dying of CaP [7]. A similar difficulty exists over the 
accurate identification of indolent disease, in which either 
radical surgery or high dose radiation treatments could be 
deferred [1,2]. 

MicroRNAs (miRNAs) hold considerable promise 
to be exploited as highly accurate and functional 
prognostic serum markers of CaP stages and drug 
responses [8-18]. They can encapsulate events within 
the tumor micro-environment and thereby overcome the 
limitations of inaccurate tumor sampling at biopsy. From 
a biostatistical perspective, given there are fewer miRNA 
than protein coding mRNA, genome-wide coverage is 
more readily achieved and avoids the statistical penalties 
typically associated mRNA genome wide testing [19]. 

Despite encouraging results from analyzing miRNA 
expression to define later stages of CaP [20, 21], it has 
yet to be demonstrated that miRNAs can be leveraged in 
combination with non-invasive clinical measures (e.g. 
PSA) at an early disease stage to identify those cancer 
patients who will rapidly recur after treatment. The 
goal of the current study was to find miRNA expression 
patterns in men with localized CaP that associate with 
the risk of early treatment failure defined by biochemical 
progression soon after RP and/or can be exploited to 
predict progression better than pathological and/or clinical 
parameters alone [22]. 

To meet this challenge we used a combination of 
cell line approaches and in silico analyses of publically 
available data derived from CaP tumors to identify 
miRNA associated with different CaP stages (Figure 1). 
The capacity of these miRNA to distinguish progression 
risks and treatment responses were subsequently measured 
in serum samples from CaP patients. Principally, we 
focused on serum collected prior to surgery, in a clinical 
cohort of men with localized CaP who underwent RP, with 
the overall goal of defining miRNA expression patterns 
associated with the risk of biochemical progression. 
Three miRNAs (miR-103, 125b and 222) were identified 
that, at the time of surgery, associated with the subsequent 
risk of biochemical progression. When combined with 
PSA levels, these miRNA generated strong predictive 
models that identified men who would experience 
biochemical progression. Critically, these models, derived 
from non- invasive measurements, were as effective 
at predicting biochemical progression as the tumor 
grade derived from the pathological examination of the 
surgically removed tumor. These findings suggest that 
serum miRNA expression patterns in combination with 
PSA levels can be exploited to stratify patients for optimal 
therapeutic response, even prior to surgery. 



2. RESULTS 

2.1. Identification of miRNA expression patterns 
in non-malignant and malignant prostate cell 
lines. 

The purpose of these experiments was to identify 
altered expression patterns of miRNAs that reflect 
different CaP stages, by agnostically analyzing miRNA 
expression levels in cell line models; specifically RWPE- 

1 vs RWPE-2 [33]; HPrlAR vs LNCaP [34]; LNCaP vs 
LNCaPC4-2 [35]. MiRNA microarray analyses of these 
three pairs of cells identified a total of 56 miRNAs as 
significantly altered (1.5 fold, /><0.05) (Figure 2A and 
Table 1). Comparison between RWPE-1 vs RWPE-2 
(encompassing CaP initiation) identified 11 significant 
changes, and of these miRNAs, 9 were up-regulated and 

2 were down-regulated in RWPE-2. Comparison between 
HPrl-AR vs LNCaP (encompassing CaP initiation & 
progression) identified 43 significant changes of which 
33 miRNAs were down-regulated and 10 miRNAs 
were up-regulated in LNCaR LNCaP vs LNCaP-C4-2 
(encompassing the ADT-CaP phenotype) revealed 16 
miRNAs that were significantly up-regulated and 3 
miRNAs were down-regulated in LNCaP-C4-2. HPrlAR 
vs LNCaP, and LNCaP vs LNCaP-C4-2 comparisons 
revealed nine common miRNAs (Figure 2A and Table 
1). Out of 11 significant changes in RWPE-1 vs RWPE-2, 
8 miRNAs were common with HPrlAR vs LNCaP and 
two with LNCaP vs LNCaP C4-2. Two miRNAs were 
differentially regulated across the three comparisons (miR- 
335 and miR-31, Figure 2 A and Table 1). Collectively, 
these findings suggest there are common and unique 
miRNA expression patterns for different stages of CaP. 

The 56 miRNAs that were identified to be 
differentially expressed by the microarray approach 
were selected for subsequent validation by Q-RT- 
PCR, using Exiqon LNA™ miRNA specific primers. 
Validated PCR primers were unavailable for 8 miRNAs, 
and 5 miRNAs failed QC and were excluded from the 
analysis. Thus, 43 miRNAs were tested for validation 
that included 10 miRNAs associated with RWPE-1 vs 
RWPE-2, 34 in HPrlAR vs LNCaP and 11 in LNCaP 
vs LNCaP C4-2 (Table 1). Q-RT-PCR vahdation of 
miRNA expression changes in the cell line pairs showed 
strong concordance with the patterns observed in the 
microarray data; 9 of 10 in RWPE-2, 33 of 34 in LNCaP 
and 4 of 1 1 in LNCaP-C4-2 showed significant changes 
similar to microarray data (>=1.5 fold,/><0.05. Table 1). 
Additionally, 17 other miRNAs in these three cell line 
combinations that failed to reach a significance level 
in microarray analysis revealed modest but significant 
differential regulation in Q-RT-PCR analysis, for example 
miR-200b in RWPE-2, miR-193b in RWPE-2 and LNCaP 
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and miR-29a in LNCaP C4-2 (Table 1). 

Subsequently, these 43 miRNAs were measured in a 
panel of 14 CaP cell models, which included the original 6 
cell lines used for the microarray analyses. 12 of these cell 



lines were derived from European American (EA) and two 
from African American (A A) patients. Relative miRNA 
expression in these cell lines was calculated compared 
to non-malignant RWPE-1 cells. The expression patterns 



1. Identification of altered miRNA in prostate cell lines 
Approach: 

Three cell lines pairs (2 isogenic) exannined in triplicate by nnicroarray 
Outcomes: 

56 significant differentially expressed nniRNA 



^ 

2. Technical validation of altered miRNA expression in cell lines 
Approach: 

Q-RT-PCR of nniRNA in 14 non-nnalignant and nnalignant cell nnodels 
Outcomes: 

Significant validation of 43 nniRNA that group cells by phenotype 





3. In Silico analyses of publically available prostate cancer data 
Approach: 

1. Connparison of cell line data with MSKCC CaP data (Taylor et al) 

2. Connparison of cell line expression differences with DNA nnethylation arrays 
from ENCODE 

Outcomes: 

1. Positive correlation of nniRNA expression in cell lines and tunnors 

2. Negative correlation between expression and CpG nnethylation 



Test prognostic capacity of serum miRNA expression 




Progressors Non-Progressors 



4. Analyses of serum miRNA expression and clinical outcome 
Approach: 

Q-RT-PCR of 43 miRNA 

Multivariate Cox proportional hazard model 

Penalized regression models (Lasso models) 

Outcome: 

3 miRNA significantly predict disease biochemical progression risk 



Figure 1: The workflow designed to identify and test the prognostic capacity of miRNA in the serum of patients with 
CaP. 
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Table 1: miRNA expression patterns identified in microarray analysis and validated by Q-RT-PCR in 
cell lines and CaP serum samples. RNA was isolated from the indicated cell lines and human tumor samples 
using Trizol and from serum samples using Qiagen miRNeasy kit with Exiqon recommended modifications. 
Microarray analysis was performed with Exiqon miRNA 5th generation expression arrays using single color 
hybridization in biological triplicates for each cell line. 384-well plates Q-RT-PCR analysis was used with 
Exiqon miRNA specific primers. All miRNAs with missing values (ct>38) were removed. Samples were then 
normalized using the overall mean miRNA expression value. MiRNA differentially expressed in microarray 
analysis and validated by Q-RT-PCR are shown in bold letters with asterisk mark. (>1.5 fold change and 
/?<0.05). ND, not detected, failed in Q-RT-PCR; NA, primers to validate not available 
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of these 43 miRNAs, detected by Q-PCR in the 14 cell 
line panel, are represented in Figure 2B. The EA cell 
lines with similar clinical features (e.g. non-malignant, 
androgen sensitive CaP and ADT-RCaP) showed related 
expression changes and, using hierarchical clustering, 
were grouped together (Figure 2B). The AA ADT-RCaP 
model (MDA PCa 2b) was comparable to the EA ADT- 
RCaP model (LNCaP C4-2). However, the AA androgen- 
sensitive model (E006AA) grouped with the metastatic 
EA cell lines (e.g. PC-3 M). These data suggest that 
altered miRNA expression magnitude and direction can 
distinguish between cell lines from different stages of CaP, 
including those derived from men of different race. 



2.2. In silica analyses of miRNA expression in 
human CaP. 

Next we tested whether differentially expressed 
miRNA patterns, identified and validated through the 
prostate cell line analyses, reflect miRNA expression in 
human CaP tumors. For this analysis we interrogated 
publically-available miRNA expression data in primary 
tumors from the MSKCC cohort [28]. Only CaP patients 
reporting as EA non-hispanics and who underwent RP or 
laproscopic RP (LRP) were used for analysis. Relative 
expression of miRNAs in primary prostate tumors (N=78) 
from EA men were calculated compared to normal prostate 
tissue from EA men (N=25). Normalized log^ expression 



A 



RWPE2/RWPE1 
(9 up, 2 down) 



B. 



7^ 




LNCaP/HPrlAR 
110 UP, 33 down) 



tkt ^ iH o m 



LNCaP-C4-2/LNCaP 
(16 up, 3 down) 




Figure 2: Identification of differentially expressed miRNAs by microarray and Q-RT-PCR in CaP cell lines. A. The 

number of differentially up and down regulated miRNAs identified in three cell line pairs representing different clinical stages of CaP. 
Non-malignant RWPE-1 and their isogenic RAS-transformed counterpart RWPE-2; non-malignant HPrl-AR and LNCaP; and LNCaP and 
their isogenic ADT recurrent variant LNCaP-C4-2 cells using Exiqon miRNA microarray. Cells were grown to mid-exponential density and 
RNA isolated prior to microarray analyses. All experiments were undertaken in triplicate. B. 43 miRNA were measured by Q-RT-PCR in a 
panel of 14 non-malignant and malignant cell lines. Expression analysis was performed using Exiqon miRNA specific primers. The relative 
expression of miRNAs in different cell lines was calculated compared to RWPEl cell line. Heat map shows log2 fold change values (green 
representing downregulation, red representing upregulation relative to RWPEl). 
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data for miRNAs were downloaded from the cBioPortal 
and compared to the RWPE-1 vs RWPE-2, and HPrlAR 
vs LNCaP microarray data. 

367 miRNA were in common between the two 
microarray platforms and their expression patterns are 
represented as scatter plots in Figure 3. These analyses 



revealed a significant correlation between expression 
patterns identified by the cell line and primary tumor 
analyses. For example, the differentially expressed 
miRNAs between RWPE-2 vs RWPE-1 and the primary 
tumor vs. adjacent normal comparison revealed a 
significant correlation of both up (e.g. miR-141) and 
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Figure 3: Scatter plot showing the correlation of miRNA expression between prostate cell lines and publically available 
CaP tumors. Log 2 fold change of miRNA expression was calculated for 78 primary CaP tumors compared to 25 matched normal tissue 
from Taylor et al study (MSKCC cohort) and compared to the log 2 fold change analyses obtained through cell line comparisons. A. 
miRNA expression primary in CaP tumors compared to RWPE2/RWPE1 cells. B. miRNA expression in primary CaP tumors compared to 
LNCaP/HPrl AR cells. Up regulated miRNAs in cell lines are shown as red circles, donwregulated miRNAs in cell lines are shown as green 
triangles. Red and green symbols show significant correlation between the data sets. Blue symbols are miRNA that were not significantly 
correlated. 
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Table 2: Patients characteristics used in case-cohort 
study of miRNA association with risk of biochemical 
progression. All serum samples were collected with 
IRB approval at RPCI. Patient serum and demographic, 
epidemiologic, clinical and pathology data were collected 
through the Data Bank and BioRepository (DBBR), 
a shared resource at RPCI. All blood specimens were 
collected from newly diagnosed men with CaP, prior to 
any treatment. All patients had at least 3 years of follow- 
up data. Biochemical progression was defined by serum 
PSA of 0.2 ng/mL or greater (obtained 6 weeks - 3 months 
postoperatively), with a second confirmatory level of PSA 
greater than 0.2 ng/mL (N=31, classified as progressors), 
all other men were classified as non-progressors (N=62). 
Non-Progressors were selected from a larger group of CaP 
patients with equivalent follow-up and frequency matched 
on age and Gleason grade sum using a bootstraping 
procedure. ^Average follow-up time for non-progressors 



Variables 


Progressors 


Non-progressors 


N 


31 


62 


Age at diagnosis 


59.2 


59.8 


PSA at Surgery 


10.7 


6.3 


Average Time to 
Progression (days)* 


624.6 


1634.6 


Gleason Grade Sum 






6 


3 


7 


7 


20 


44 


8 


2 


2 


9 


6 


9 



down (e.g. miR-31) regulated miRNA (p<0.0036). 
Similar significant findings (p<0.0001) were identified 
by comparison of LNCaP cells compared to HPrlAR 
and miRNA changes in primary tumors (e.g. miR- 
222). This in-silico validation suggested that miRNA 
expression patterns identified through cell line analysis are 
comparable to expression in human CaP tumors (Figure 
3). 

2.3. Serum miRNA levels at the time of RP 
significantly associate with the risk of biochemical 
progression. 

The ability of the 43 miRNAs, identified through 
cell line and in silico analyses, to predict the risk 
biochemical progression was tested in a well-defined 
clinical cohort of 99 men with primary CaP, who were 
frequency matched on age and Gleason grade sum. 
Six serum samples failed Q-PCR quality control and 



Table 3: Hazard of biochemical progression. 

Multivariate Cox proportional hazard models were 
used to determine the miRNA, pathological and clinical 
features that were significantly associated with hazard of 
biochemical recurrence. 



Model 


HR 


Lower 
95% CI 


Upper 
95% CI 


Model 
p-value 


PSA 


1.05 


1.02 


1.08 


0.001 


Gleason Grade Sum 


1.22 


0.82 


1.80 


0.32 




\j .yy 


\j .y D 
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BMI 
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1.12 
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hsa.miR.222 


2.80 


1.29 


6.10 


0.009 


hsa.miR.103 


0.41 


0.21 


0.79 


0.008 


hsa.miR.125b 


1.79 


1.10 


2.91 


0.018 


hsa.miR.222, 
adjusted for PSA 


2.60 


1.17 


5.76 


0.0004 


hsa.miR.103, 
adjusted for PSA 


0.41 


0.20 


0.82 


0.0002 


hsa.miR.125b, 
adjusted for PSA 


1.78 


1.08 


2.95 


0.0006 



were removed from subsequent analysis; the patient 
characteristics of these 93 samples are described in Table 
2 (62 patients who did not progress during follow-up and 
31 patients who progressed). 9 miRNAs were undetected 
in serum (ct > 38) and were removed from the further 
analysis. Of the 34 miRNAs remaining that were serum 
detectable, 3 were found to be differentially expressed in 
the serum of men who progressed compared to those who 
did not (miR-103, miR-125b, mir-222). The proportion 
of men who progressed is shown for these three miRNAs 
serum expression values and PSA quartiles in Figure 4A. 
The proportion of progressors is significantly different 
between first and fourth miRNA quartile (/7<0.0001) and 
with PSA level (p<0.001). 

We then measured the hazard of biochemical 
progression attributable to clinical and pathological 
characteristics, as well as serum expression for each of 
these three miRNA (normalized ct values) (Table 3). We 
found a modest but significant association with PSA such 
that for each unit increase in PSA at diagnosis there was a 
5% increase in hazard of biochemical progression. Given 
the frequency-matching schema, in which patients in the 
two groups were matched based on Gleason grade sum 
and age, it is not surprising that Gleason grade sum was 
insignificant as was age and BMI at time of diagnosis. The 
hazard associated with each unit increase in normalized ct 
value of miR-222 was 2.8 (1.3 - 6.1 95% CI) and for miR- 
125b it was 1.8 (1.1 - 2.9 95% CI). In contrast, each unit 
decrease in miR-103 associated with an almost 2.5 fold 
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decrease in hazard of biochemical progression (HR=0.41 
with 0.2-0.8 95% CI, Table 3). 

To make the clinical data more suitable for analyses, 
we separated it into discrete categories (quartiles). Again, 
not surprisingly, quartile analyses of these miRNAs (first 
and fourth) remain significant at a similar magnitude. 
Quartile four values for miR-222 and miR-125b 
revealed a significant increase in hazard of biochemical 
progression compared to quartile one; CaP patients with 
miR-103 values in quartile one showed a decrease risk 
of progression compared to men in quartile four (Figure 
4B). These three miRNA were highly correlated with one 
another: miR-222 and miR-125b positively correlated {p 
=0.74); miR-222 and miR-103 negatively correlated {p = 
-0.37), as were miR-125b and miR-103 (p = -0.56). 

Importantly, the miRNA values and clinical 
variables represent different contributions to hazard of 
biochemical progression, as reflected by the fact that no 
significant correlations were found between PSA, Gleason 
score, age, BMI and any of the three significant miRNA. 
The best fitting hazard model included both miR-103 
and PSA (p<0.0002), although all three miRNAs when 
combined with PSA were of similar significance levels 
(Table 3). 

A. 



We also measured serum expression of the 34 
miRNA used above in two small cohorts of men who had 
progressed after RP and were undergoing treatments for 
non-localized metastatic stage. Specifically, these were 
men with hormone sensitive prostate cancer (HSPC, 
N=17) treated with anti-androgen therapy, and men with 
ADT-RCaP being treated with chemotherapy (N=14). 
In the HSPC group, 11 miRNA differed between time 
points, the most significant, miR-320a (log^ fold change 
=-0.92, p<.0001. Supplementary Table 1). Of the miRNAs 
identified as differentially expressed in men with localized 
disease, miR-222 was also differentially expressed 
between the time points (log^ fold change=-0.37,/><.04). 
In the ADT-RCaP group, we measured association of 
miRNA with hazard of death (5/14 patients died). While 
none of these 3 miRNA were significant alone, a model 
including PSA and miR-103 was suggestive of statistical 
significance (p<0.06), although due to the small sample 
size the hazard ratio confidence intervals were wide 
(HR=2.65 with a 95% CI of 0.55-15.5). These findings 
also added confidence to the identification of these 
miRNA in the RP cohort as being associated with disease 
progression. 
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Figure 4: A. Distribution of miRNA serum expression values within the RPCI cohort MiRNA expression levels were divided into quartiles 
for the entire cohort. For each quartile the proportion of the total number of samples (y axis) is shown for progressors and non-progressors. 
B. Time to biochemical progression analyses for miRNA. Quartile analyses of these miRNA (first and fourth) remain significant at a 
similar magnitude, with quartile four values for miR-222 and miR-125b showing significant increase in hazard of biochemical progression 
compared to quartile one; CaP patients with miR-103 values in quartile one showed a decrease risk of progression compared to men in 
quartile four. Time is measured in days from surgery. 
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2.4. Correlation of miRNA expression in tumor 
and serum. 

Finally, we also measured expression of the 34 
miRNAs in matched serum and primary tumor from 
a limited subset of ten CaP patients from RPCI whose 
serum samples were analyzed as part of the biochemical 
progression analyses. All three miRNA, miR-103, mir- 
125b and miR-222, showed relatively similar expression 
levels in the tumor and the serum in each of the ten 
patients suggesting that, for these miRNA at least, the 
levels detected in serum accurately reflect the expression 
of miRNA in their tumor counterparts. This was not 
universal for all miRNA. 

2.5. Generating predictive models of biochemical 
progression. 

To assess the predictive value of the miRNA 
for biochemical progression, we undertook penalized 
regression analyses with leave one out cross validation 
with all 34 serum detected miRNAs between the two 
groups, PSA, age and pathological grade [31, 32]. The 
AUC for the restricted lasso model containing PSA, age 
and pathological grade was 0.55 (95% CI, 0.38-0.64, blue 
line in Figure 6). The unrestricted lasso model again, 
considering all miRNA and clinical parameters, contained 
covariates PSA, and three miRNAs out of 34 tested (miR- 
103, miR-125b, and miR-222) with AUC 0.64 (95% CI, 
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0.55-0.79, Red line in Figure 6). 

This finding underscores the discovery of significant 
changes in expression of these three miRNA between men 
who progressed compared to those with stable disease. 
Although there was a difference between the predictive 
value of the AUC, a bootstrap test of significance 
comparing restricted versus unrestricted lasso model 
indicated it was not significant at the 0.05 level (p=0.25). 
Similar AUC estimates were obtained using two separate 
approaches [31, 32]. Together, these analyses revealed that 
serum expression levels of three miRNAs in conjunction 
with PSA are equally able to predict men who are at 
risk of experiencing treatment failure as the analyses 
of pathological parameters obtained after the radical 
procedure of surgery. 

3. DISUSSION 

There is considerable enthusiasm to exploit serum 
miRNA expression in diagnosis and prognosis in cancer. In 
CaP, miRNA patterns have been associated with extremes 
of disease progression, for example comparing either 
tumor or metastatic disease [36, 37] to non-malignant 
conditions [18, 28, 38-40]; or comparing prostate cancer at 
the time of surgery to benign prostatic hyperplasia [37, 41, 
42]; or comparing different tumor grades [43]. Although 
there is a rapidly expanding literature on CaP miRNA 
serum expression patterns (Supplementary Table 2), the 
current study is significant for several reasons. Firstly, the 
current approach is the first study to establish that miRNA 
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Figure 5: miRNA expression in matched tumor and serum samples. Expression of miRNA was analyzed in paired serum and 
tumor samples from 10 CaP patients. miRNA expression was performed using Exiqon miRNA specific assay by Q-RT-PCR. Raw ct values 
were mean centric normalized then difference between serum and tumor samples were calculated (dCT^ -dCT^ ). 
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in combination with non-invasive clinical parameters 
could be used to predict early CaP progression without 
knowledge of the underlying tumor grade. Secondly, the 
use of serum miRNA is as effective, and suggestive of 
being better, than analyses of the pathological parameters 
alone that are recovered following surgery. 

We began with a thorough analyses of cell line 
models, reasoning that cell lines are more pure cell 
populations compared to primary tissue, and therefore 
would yield more accurate differential miRNA expression 
patterns to carry forward to a human cohort. Microarray 
and Q-RT-PCR approaches in a 14 prostate cell line panel 
supported the concept that the expression patterns of 43 
miRNA could accurately distinguish of non-malignant, 
androgen sensitive and ADT-RCaP cell types. This 
miRNA panel includes miRNAs previously associated 
with resistance to chemotherapies (miR-130a) [44] and 
poor progression free survival (miR-200c and miR-141) 
in CaP [45]. These also included miRNAs associated 
with regulation of proliferation, metastasis and poor 
progression free survival [46-48]. Pathway enrichment 
analyses [29] of these differentially expressed miRNA 
identified significant targeting of TGF-P and WNT 
networks (reviewed in[49, 50]) (Supplementary Table 
3). For example, the comparison of RWPE-2 to RWPE- 
1 identified 3 miRNAs (e.g. miR138, 200c and let-7i) 



Lasso, restricted vs. unrestricted models 




1.0 0.8 0.6 0.4 0.2 0.0 



Specificity 

Figure 6: Predictive models of biochemical progression. 

To assess predictive value of the significant miRNAs, penalized 
regression analyses was performed including all miRNA, 
PSA, age and grade. The AUC for the restricted lasso model 
containing PSA, age and grade was 0.52 (95% CI 0.38-0.64, 
Blue line). The unrestricted lasso model contained covariates 
PSA, miR-103, miR-125b, and miR-222 with AUC 0.67(95% 
CI, 0.55-0.79, Red line). 



that targeted thirty-four genes on the TGF-P pathway, 
including TGBRl and SMAD4 [51, 52]. Interestingly, 
loss of Smad4, a downstream regulator of TGF-P, drives 
invasive CaP in co-operation with Pten depletion [53]. The 
cell lines represent different stages of disease progression 
and it is interesting to note that the miRNA expression 
patterns continue to change between the different cell 
types and disease states; miR-138 was up-regulated in 
RWPE-2 compared to RWPE-1, but down-regulated in 
the androgen sensitive cancer cell lines (Table 1). This 
re-enforces the concept that unique and shared miRNA 
expression patterns (both in terms of magnitude and 
direction) may precisely reflect CaP disease states, and 
determine the risk of CaP progression. 

Interestingly, epigenetic events are implicated in 
the control of miRNA expression [54], (reviewed in[55]). 
To examine the basis for the altered miRNA expression 
we considered the impact of DNA CpG methylation. 
We mined ENCODE data sets of genome wide DNA 
methylation patterns in LNCaP compared to PrEC non- 
malignant prostate cells to examine the distribution of 
methylation at CpG sites within ±5kb of the 35 miRNA 
found in the current study to be altered in LNCaP as 
compared to RWPE-1 cells. 23 miRNA displayed an 
inverse relationship between DNA CpG methylation 
status and miRNA expression levels, suggesting that 
CpG methylation impacts the expression of miRNA 
in CaP to a significant extent (Supplementary Figure 1, 
Supplementary Table 4) and may be associated with tumor 
initiation. 

Of the 43 miRNA identified in the cell line 
comparison only 34 were adequately detected in the 
serum of CaP patients, and of these 3 miRNAs (miR- 
103, miR-125b, miR-222) were differentially expressed 
at the time of surgery between patients who were 
otherwise indistinguishable in terms of age, PSA levels, 
tumor stage derived from biopsy material but diverged in 
terms of displaying either stable disease or biochemical 
progression. Models were constructed to measure the 
relationship between serum miRNA expression, the hazard 
of progression and the ability of the miRNA to predict 
progression alone and in conjunction with PSA, stage, 
grade and BMI. Irrespective of the modeling approach, 
our results consistently identified that differential 
expression of miR-103, miR-125b and miR-222 predicted 
biochemical progression. Of these, miR-222 [14, 56-61] 
and to a lesser extent miR-125b [21, 41, 42, 59, 62] appear 
to play an important role in prostate biology, whereas to- 
date miR-103 has not been associated in CaP, although has 
been associated with other cancers [63-65]. 

MiR-222 is transcribed as a miR cluster with miR- 
221 from chromosome Xllp.3. Felli et al. (2005) first 
described miR221/-222 as integral regulators of in vitro 
erythropoiesis[66] and deregulated expression of the 
cluster is implicated in several cancers [67-75] suggesting 
that these miRs have oncogenic potential. In general. 
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oncogenic properties of miR-221 and miR-222 are 
attributed to their control of the cychn dependent kinase 
(Cdk) inhibitors p27KIPl and p57KIP2, and thus control 
of Gl to S phase transition[56, 57] PI3K and PTEN 
signaling [14] and other targets implicated in malignant 
transformation, including CX43[68], RECK[71], and 
ERa[76, 77]. 

Tissue miR-222 expression was also recently found 
to predict progression free and overall survival in gastric 
carcinoma patients [73], and miR-222 detection in urine 
was shown to detect bladder cancer with high accuracy, 
thus providing a potential noninvasive diagnostic 
biomarker in this setting, as well as predicting for 
recurrence, progression, and overall survival[75]. In CaP, 
Lin et al. (2011) observed that tumors (Gleason grade > 
7) had elevated expression of miR-221 and miR-222 
compared to less aggressive tumor tissues (Gleason grade 
< 7), further suggesting an oncogenic role of these miRs in 
CaP development[58]. Interestingly, an assessment of 40 
early stage (T2a/b) prostatectomy specimens revealed loss 
of miR221/-222 expression in microdissected malignant 
tissue compared to non-involved tissue controls [5 9], 
perhaps indicating that expression levels are differentially 
modulated at divergent stages of disease and evidence also 
suggests that miR221/-222 is involved in ADT-RCaP[60]. 
Most recently, a miRnome wide scan found miR-222 to be 
the most significantly differentially expressed miR when 
comparing 20 matched pairs of microdissected tissue 
samples of prostate cancer and non-tumor tissue[78]. 

MiR- 125b was also elevated in men with 
progressive disease in the current study although the 
function is ill-defined as it has been reported to function 
as oncogene or tumor- suppressor gene in different cancer 
types or cell lines [79]. One of the earliest studies on miR 
in CaP identified six miR including miR- 125b and miR- 
143 to be upregulated in metastatic CaP serum samples 
as compared to normal individuals [17]. Expression 
of miR- 125b in serum of CaP patients is reported to be 
upregulated as compared to normal controls[21] whereas 
other studies reported it to be downregulated in CaP as 
compared to normal or BPH samples[41, 42, 59]. MiR- 
125b regulates cell proliferation in prostate cancer cell 
lines cells [62], and it has suggested to be upregulated by 
androgen signaling[80]. Functionally in CaP, miR- 125b 
has been reported to target BAK1[^0] (a pro-apoptotic 
member of the BCL-2 gene family) and EIF4EBP1 [41] 
(Eukaryotic translation initiation factor 4E-binding 
protein 1 , a gene that encodes one member of a family of 
translation repressors proteins) and the transcriptional co- 
repressors NC0R2/SMRT [81]. 

Our finding of reduced miR- 103 in the serum of 
men whose disease progressed and association with 
hazard of progression is a novel one in CaP. The aberrant 
overexpression of miR- 103 has been reported in many 
cancers including endometrial[63], bladder[64] and 
breast cancers[65]. Increased expression of miR- 103 was 



associated with metastasis and poor outcome in breast 
cancer[65]. Functionally miR- 103 targets Dicer and 
thus can attenuate miR biosynthesis. MiR- 103 has also 
been shown to regulate cancer metastasis by regulating 
EMT and inhibition of miR- 103 opposes migration and 
metastasis[65]. 

Tumor-serum correlation analysis of miRNA 
expression in a subset of these patients revealed variations 
in detectible expression suggesting parallel processes 
which govern serum miRNA detection [82-84] including 
active miRNA export where the serum levels were high, 
and passive release for miRNA with lower and poor 
correlations between serum and tumor expression levels. 
For the 9 miRNA that were undetectable in serum, despite 
robust expression in tumors, it suggests that there are 
differential secretion of miRNA, or selective degradation 
in the serum, or that secreted miRNA profile reflects a 
more complex interaction of different cell types within the 
tumor micro-environment. 

The potential for altered tumor and serum 
relationships was also highlighted in our preliminary 
validation studies. In this case we examined the 
association of miR- 103, miR- 125b and miR-222 
with biochemical progression in the 78 white non- 
hispanic patients from the MSKCC data set. Patients 
were separated into those with high or low miRNA 
expression using a median expression level threshold, and 
proportions of each subgroup to experience biochemical 
progression over time were calculated as Kaplan-Meier 
curves (Supplementary Figure 2). MiR- 125b expression 
associated with time to biochemical progression (p<0.06) 
in this cohort. Interestingly, low CaP tumor expression 
of miR- 125b predicts increased overall probability of 
progression, and is opposite of what we observed in 
the serum of CaP patient, where elevated expression 
significantly predicts increased probability of biochemical 
progression. However, active secretion of miR- 125b has 
been reported in serum upon castration in mice and in 
extracellular media upon androgen blockade in LNCaP 
cells with simultaneous loss of expression in cellular 
compartments [81], most likely downstream of AR 
actions. Therefore these trends may indicate that patients 
with the most aggressive disease have most active export 
from the tumor (and therefore lowest tumor levels) and 
the highest levels of miR- 125b in the serum as a result. 
Analysis of the expression of miR- 125b in the matched 
tumor serum samples was supportive of this relationship 
also. The expression levels in the matched tumor and 
serum samples was suggestive of lower tumor expression 
and higher serum levels in the progressors compared to the 
non-progressors (Supplementary Figure 2). 

Recently, Selth et al., [13] demonstrated there may 
be evidence for association of serum expression of miR- 
146-3p and miR- 194 with reduced time to biochemical 
progression. Specifically, the authors suggested that 
high levels of miR-146b-3p expression, pre-operative 
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PSA and seminal vesicle invasion, obtained following 
invasive surgery, were predictors of a reduced recurrence 
free interval. Furthermore, this miRNA was either very 
low or not detectable in greater than one half the samples 
in either of the two publically available CaP tumor 
cohorts examined (MSKCC or Erasmus) precluding 
further evaluation. When considering only miR-194 
expression levels values, quartile one versus four, there 
was a significant association with time to biochemical 
progression in the MSKCC and the Selth et al. cohort, 
respectively. 

However, miR-194 expression did not show 
significant association with time to biochemical 
progression in the Erasmus cohort. Interestingly, in the 
MSKCC cohort the CaP patients with metastases, who 
comprise a majority of the low miR-194 group, drove 
the statistical differences for these comparisons. The 
prominent role of metastatic cancer (versus primary) 
implies that miR-194 could be a marker of other events 
such as general inflammation rather than specifically 
biochemical progression. Furthermore, the unique 
contribution of the miRNA expression to biochemical 
progression was hard to establish, in that both miR-194 and 
miR-146-3p were weakly correlated with PSA and highly 
correlated with age and surgical margins, respectively; the 
later association perhaps explaining the uniqueness of the 
authors miR-146-3p finding. Interestingly neither miRNA 
was significantly differentially expressed in our cell line 
data, and thus was not investigated in the current human 
CaP cohort. 

This is neither surprising nor uncommon. For 
example, diagnostic gene expression tests that predict 
breast cancer progression risks include different gene 
classifiers but are equally effective and clinically approved 
[85, 86]. These considerations aside, this recent study 
[13] combined with the data in the current study provide 
evidence that indeed miRNA expression levels in serum 
can be found to associate with biochemical progression 
and that this association reflects critical events in the 
tumor. 

There is intense interest in the exploitation of 
serum miRNA patterns as prognostic markers to predict 
progression risk directly. Uniquely, the present study 
revealed that miRNA expression patterns in serum, 
collected before surgery, and combined with PSA are 
better able to predict early progression risk than merely 
clinical variables alone (Table 3). That is, the current 
study remains (to the best of our knowledge) the first to 
test miRNA prognostic capacity using a cohort of matched 
patient serum samples. In this capacity we propose that 
the current study will encourage the development of 
other prognostic miRNA signatures, either in the same 
or in different prostate cancer disease settings. The 
generalizable nature of individual miRNA to predict risk 
will only emerge from the meta-analyses of such studies. 

In turn, these findings have significant potential 



to reduce the impact of co-morbidities associated with 
surgery [87, 88] and to identify patients for whom 
aggressive therapies are warranted. The prospect of 
being able to exploit serum miRNA expression, as a PCR 
based assay, either alone or combined with other serum 
markers (e.g. PSA) has a profound capacity to improve 
the clinical management of prostate cancer by widespread 
application of a cost-effective, robust and non-invasive 
serum measurement. 

4. MATERIALS AND METHODS 



4.1. Biological analyses: 

4.1.1. Cell Lines: 

Prostate cell lines included those derived from 
European American men; normal prostate epithelial 
cells (RWPE-1, HPrl, HPrlAR and PNT2), early stage 
CaP cells (RWPE-2), androgen sensitive CaP (LNCaP), 
Androgen deprivation therapy (ADT) recurrent CaP (ADT- 
RCaP) (LNCaP-C4-2), metastatic CaP (DU145, PCS and 
PC3M) and CaP xenograft cell lines (22Rvl and LAPC4). 
Cell lines from African American men; clinically localized 
CaP (E006AA)[23] and ADT recurrent CaP (MDA PCa 
2b) [24]. All cells were authenticated at Roswell Park 
Cancer Institute (RPCI) immediately prior to the start of 
the study. 

2.1.2. RPCI Patient Samples: All serum samples 
were collected with IRB approval at RPCI. Patient 
serum and demographic, epidemiologic, clinical and 
pathology data were collected through the Data Bank 
and BioRepository (DBBR), a shared resource at 
RPCI. All blood specimens were collected from newly 
diagnosed men with CaP, prior to any treatment, and were 
immediately processed and stored in liquid nitrogen within 
one hour of blood draw [25]. All men had newly diagnosed 
clinically localized CaP and had no prior history of either 
another cancer or other treatment for CaP at the time of 
blood collection. Standard prognostic variables included 
clinical (cTNM) and pathological (pTNM) stage, Gleason 
score, and PSA. All patients had at least 3 years of follow- 
up data. Biochemical progression [5] was defined by 
serum PSA of 0.2 ng/mL or greater (obtained 6 weeks - 
3 months postoperatively), with a second confirmatory 
level of PSA greater than 0.2 ng/mL (N=33, classified 
as progressors), all other men were classified as non- 
progressors (N=66). The progressors and non-progressors 
were selected from a larger group of CaP patients who 
attended urology clinics at RPCI, and for whom there 
was detailed follow-up (N>450) extending over an 8 year 
period. From these samples the test cohort was selected 
by being frequency matched on age and Gleason grade 
sum using a bootstrapping procedure. For 10 of these 
selected serum samples, matched tumor samples were also 
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available. 

4.1.3. Microarray and Q-RT-PCR analysis: 

RNA was isolated from cell lines and human 
tumor samples using Trizol and from serum samples 
using Qiagen miRNeasy kit with Exiqon recommended 
modifications. Microarray analysis was performed with 
Exiqon miRNA 5th generation expression arrays using 
single color hybridization in biological triplicates for 
each cell line. Q-RT-PCR analysis was performed using 
Exiqon miRNA specific primers and reagents in 384 well 
plates as recommended by Exiqon. All miRNAs with 
missing values (ct>38) were removed. Samples were then 
normalized using the overall mean miRNA expression 
value [26]. 

4.2. Statistical analyses: 

4.2.1. Microarray: 

Normalized expression of miRNAs was compared 
across pairs of cell lines using linear mixed effects models 
as implemented in LIMMA[27]. Differentially expressed 
miRNAs were identified using p-value <0.05 and fold 
change >1.5 cut-off. 

4.2.2. Q-RT-PCR: 

Differentially expressed miRNAs compared to 
RWPE-1 were normalized using 5S rRNA (Exiqon 
reference gene panel) and measured with two-tailed t-test 
(p-value <0.05 and fold change >1.5 cut-off). Hierarchical 
clustering was used to assess if selected miRNAs could 
be used to differentiate normal epithelial, CaP progression 
and metastatic cell lines. 

4.2.3. In silica validation: 

Publically available data was used to help determine 
if miRNA expression identified from comparing pairs of 
cell lines (RWPE-1 vs RWPE-2; HPrl-AR vs LNCaP) 
correlated with miRNA expression between tumor and 
normal samples in a subcohort of localized tumors from 
the Memorial Sloan-Kettering Cancer Center (MSKCC) 
archive [28]. Specifically, 78 self-reporting White non- 
Hispanic men with confirmed locally-confined CaP tumors 
(with 22 corresponding matched normal tissue) obtained 
from radical prostatectomy were selected for further 
analyses [28]. Information is publically available on a total 
of 368 miRNAs of which 367 were common between the 
microarrays analyzed at RPCI and MSKCC. 

4.2.4. Network analyses: 

MiRNA-targeted networks were identified by 
pathway enrichment analyses using Diana-mirPath[29]. 

4.2.5. miRNA serum tumor correlation: 

Q-RT-PCR was undertaken in 10 CaP matched 
serum-tumor pairs and Pearson's correlation calculated. 



4.2.6. RPCI Patient Samples: 

The distribution of clinicopathological 
characteristics was compared among the groups using 
analysis of variance (ANOVA) for continuous variables 
and the x^-test for categorical variables. Continuous 
variables tested included age at diagnosis, preoperative 
PSA and body mass index (BMI); Gleason grade 
sum (greater than 7 vs less than 7) was examined as a 
categorical variable. 

Multivariate Cox proportional hazard models were 
used to determine the miRNA, pathological and clinical 
features that were significantly associated with hazard of 
biochemical progression [30]. 

Penalized regression models (Lasso models using 
the R package glmnet [31, 32]) were used to determine 
what combinations of all measured miRNA, clinical 
and pathological variables best predicted biochemical 
progression. These models were also used to determine 
if either individual or joint combinations of miRNA were 
at least equivalent and/or better predictors of biochemical 
progression than clinical variables alone. Specifically, an 
unrestricted model allowing any parameter to enter/exit 
the model was compared to a restricted model that forced 
PSA, Gleason grade sum, and age to remain in the model. 
The lambda parameter for each model was estimated using 
10-fold cross-validation and binomial deviance distance 
measure. 
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